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Resonances of the time evolution (Frobenius-Perron) op- 
erator V for phase space densities have recently been shown 
to play a key role for the interrelations of classical, semiclas- 
sical and quantum dynamics. Efficient methods to determine 
resonances are thus in demand, in particular for Hamiltonian 
systems displaying a mix of chaotic and regular behavior. We 
present a powerful method based on truncating P to a finite 
matrix which not only allows to identify resonances but also 
the associated phase space structures. It is demonstrated to 
work well for a prototypical dynamical system. 

PACS numbers: 05.45.-a, 05.45.Mt, 05.20.-y 



Effectively irreversible behavior of classical Hamilto- 
nian systems can be elucidated by studying the phase 
space density and its propagator, the Frobenius-Perron 
operator V. Due to Liouville's theorem V can be repre- 
sented by an infinite unitary matrix whose spectrum lies 
on the unit circle in the complex plane. Nevertheless, 
means and correlation functions of observables can relax 
(see figure |^c) with damping factors known as (Ruelle- 
PoUicott) resonances oiV. These resonances have 
recently attracted attention, e.g. in a superanalytic ap- 
proach to universal fluctuations in quantum (quasi-) en- 
ergy spectra which originated from the physics of disor- 
dered systems. In that approach the Frobenius-Perron 
resonances constitute a link between classical and quan- 
tum chaos . There is even a recent experiment where 
quantum fingerprints of Ruelle-PoUicott resonances are 
identified [Q. To further clarify the interrelations be- 
tween classical, semiclassical, and quantum behavior, a 
practical scheme to actually determine classical reso- 
nances is called for which is free of restrictions of pre- 
vious investigations, like hyperbolicity, one-dimensional 
(quasi-) phase space or isolation of the phase-space re- 
gions causing intermittency. 

Our thus motivated quantitative investigations into 
Hamiltonian systems with mixed phase space, a still 
largely unexplored area of great interest and promise, 
lead us to the discrete unimodular Frobenius-Perron 
eigenvalues (with eigenfunctions localized in islands of 
regular motion around elliptic periodic orbits, see figure 
|l|a) and to resonances, smaller than unity in modulus 
(with eigenfunctions localized on the unstable manifolds 
of hyperbolic periodic orbits, see figure §a-c). Both the 
discrete spectrum and the resonances are determined by 
diagonalizing truncated Frobenius-Perron matrices 7^(^) 



and studying the cut-off dependence of its eigenvalues. 
Using the information about eigenfunctions, we then 
reproduce resonances by the so-called cycle expansion 
of periodic-orbit theory and furthermore through decay 
rates of correlation functions. We should add that sim- 
ilarly motivated but technically different (not involving 
eigenfunctions and employing external noise) efforts to 
determine resonances can be found in Refs. |]8|-|l0|. 

As a prototypical dynamical system we have employed 
the kicked top, i.e. a periodically kicked angular momen- 
tum J = [j svaO cos Lp, j sm. 6 sin ip, j cos 9) of conserved 
length j whose phase space is the sphere J^/j^ = 1; we 
confront a single degree of freedom with the "azimutal" 
angle </> as the coordinate and the cosine of the "polar" 
angle the conjugate momentum. The dynamics is spec- 
ified as a stroboscopic area preserving map M on phase 
space. It consists of rotations Rz{Pz), RyiPy) about the 
y— and z— axes and a "torsion" , i.e. a nonlinear rotation 
Rz {t cos 6) about the z— axis which changes (phy t cos 6, 



M^Rz{t cos e)Rz{l3z)Ry{l3y). 



(1) 



The equivalent map of the phase space density p is 
brought about by the Frobenius-Perron operator V, 



Vpicose, ip) = p{Ar'^{cose,(p)) 



(2) 



We keep f3z ~ f3y — 1 fixed and vary the torsion constant 
T, starting with the integrable case t — 0. Increasing 
values of r bring about more and more chaos until for 
r > 10 elliptic islands have become so small that they 
are difficult to detect. We shall focus on t = 4 (roughly 
90% of the phase space dominated by chaos) and r = 10 
(more than 99% chaos). 

A Hilbert space of phase space functions on the sphere 
is spanned by the spherical harmonics Yim{0,ip) with 
I = 0,1,2,... and |m| < I. These functions are or- 
dered with respect to phase space resolution by the in- 
dex I: if all Yim with < Z < Imax are admitted phase 
space structures of area roughly 4:Tr/{lmax + 1)^ can be 
resolved. If we so truncate the infinite Frobenius-Perron 
matrix Vim,i'm', we (i) destroy unitarity, (ii) restrict the 
spectrum to A'' = {Imax + 1)^ discrete eigenvalues whose 
moduli cannot exceed unity, and (iii) renounce the reso- 
lution of phase space structures of linear dimension be- 
low / {Imax + !)• Upon diagonalizing the truncated 
A^ X A^ matrix T't^) and increasing N we find the "newly 
born" eigenvalues close to the origin while the "older" 
ones move about in the complex plane. "Very old" ones 
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eventually settle for good. If the classical dynamics is 
integrable (r = or /3j, = 0), the asymptotic large— iV 
loci are back to the unit circle, where the full V has its 
spectrum. But not so for a mixed phase space: while 
some eigenvalues of ■p(^) "freeze" with unit moduli, oth- 
ers come to rest inside the unit circle as iV — > cx3. Table 
H illustrates how non-unimodular eigenvalues found for 
the kicked top with r = 10 at Imax — 40 remain in their 
positions as Imax is increased to Imax — 50, 60 and 70. 

We could pass over such findings and speak of the dan- 
ger of tampering with infinity, were there not good rea- 
sons for and a physical interpretation of the existence of 
such stable non-unimodular eigenvalues. 

The following qualitative argument suggests the per- 
sistence of non-unimodular eigenvalues as — > oo for 
non-integrable dynamics. In contrast to regular motion, 
chaos brings about a hierarchy of phase space structures 
which extends without end to ever finer scales. A trun- 
cated Frobenius-Perron operator T'^^) must reflect the 
flow of probability towards the unresolved scales as a loss, 
however large the cut-off N may be chosen. 

Arguments from perturbation theory ||l^,|l^ indicate 
that any non-unitary approximation to a unitary opera- 
tor with continuous spectrum has some eigenvalues in po- 
sitions near (non-unimodular) resonances of the unitary 
operator, i.e. poles of the resolvent in a higher Rieman- 
nian sheet. The perturbation series for such an eigenvalue 
does not converge but produces, with increasing order, a 
sequence of points concentrated in the neighborhood of 
the respective resonance. It is intuitive to interpret the 
freezing of non-unimodular eigenvalues (which need not 
be a strict convergence) as analogous to the "spectral 
concentration" known from perturbation theory. 

To find further evidence for our interpretation of frozen 
eigenvalues as resonances we have looked at the eigen- 
functions of p(^), with the following salient results. 
Eigenvalues freezing with unit moduli have eigenfunc- 
tions located on elliptic islands of regular motion sur- 
rounding elliptic periodic orbits in phase space. Such 
islands are bounded by invariant tori which form im- 
penetrable barriers in phase space. We can thus expect 
the function constant inside the elliptic islands around 
a p— periodic orbit and zero outside to be an eigenfunc- 
tion of V with eigenvalue unity. Similarly reasoning we 
expect, for p > 1, the pth roots of unity to arise as eigen- 
values as well; their eigenfunctions should have constant 
moduli and be invariant under . For the kicked top 
with r = 4 and Imax — 60 an eigenfunction with the 
eigenvalue 0.9993, i. e. almost at unity, is shown in figure 
|l|a. It is localized on the three islands around an ellip- 
tic orbit of period three (see figure |^b) and does have 
the two expected partners. We have indeed found frozen 
eigenvalues near the p-th roots of unity and their eigen- 
functions localized near elliptic period-p orbits for p up to 
6; without much further effort such signatures of higher 
periods could be identified. 



Now on to the eigenvalues freezing with moduli smaller 
than unity. Once such freezing has been observed the cor- 
responding eigenfunction has approached its final shape 
on the resolved phase space scales. The eigenfunctions 
are sharply localized around unstable manifolds of hy- 
perbolic periodic orbits, ones with low periods at first 
since these are easiest to resolve; but with growing Imax 
more complex orbits of higher periods appear in the "sup- 
port" of eigenfunctions. Even though all periodic orbits 
contributing to the structure of an eigenfunction have 
similar stability coefficients and even though the latter 
do describe the rate of mutual departure of neighboring 
trajectories it would be too naive to simply identify reso- 
nances with stability coefficients; we shall rather have to 
resort to cycle expansions further below. 

Just as for the eigenvalues there is no strict conver- 
gence of the eigenfunctions. With increasing resolution 
new structures on finer scales become visible, in corre- 
spondence with the infinitely convoluted shape of the 
unstable manifolds (see figure §a-b). Since no finite 
approximation T't^) accounts for arbitrarily fine struc- 
tures one encounters the aforementioned loss of proba- 
bility from resolved to unresolved scales. Not even in the 
limit N ^ oo can the unitarity of V be restored: rather, 
the eigenfunctions tend to singular objects outside the 
Hilbert space, in tune with a continuous spectrum of V . 

The reader may have noticed that all eigenvalues in 
table I are real or almost imaginary. In fact, all eigen- 
values we have identified as frozen inside the unit circle 
have phases corresponding to those of roots of unity, a 
fact demanding explanation. Clearly, since r^^'> is real 
the eigenvalues are either real or come in complex conju- 
gate pairs, but no other phase than zero is distinguished 
by that argument. Again, the eigenfunctions offer further 
clues. We find that the phases of the complex eigenvalues 
are determined by the length p of the shortest periodic 
orbit present in an eigenfunction / as those of the p-th 
roots of unity. The following intuitive argument indicates 
that this is to be expected. 

Assume an eigenfunction is mostly concentrated 
around a shortest unstable orbit with period p as well 
as a longer one with period p' . Denote by (5p_„ a "charac- 
teristic function" which is constant near the n-th point of 
the period-p orbit, n = 1 . . .p, and zero elsewhere. The 
truncated Frobenius-Perron operator 'P'^-' maps 5p „ into 
1 with the real positive factor rp 
smaller than unity accounting for losses, in particular 
to unresolved scales. Independent linear combinations of 
the Sp^n can be formed as fpk = e'^'^'^"^^'^p,ri with 

k = 1 . . .p. Now consider a sum of two such functions, 
9 = fpk + fp'k' , and apply T'^^) . For g to qualify as an ap- 
proximate eigenfunction we must obviously have Vp « rpi 
and k/p = k'/p'. But then indeed V'^^^'g « Vpe^'^'^'^IP g 
and [P^^^Yg « r^g. The phase is thus dictated by the 
shortest orbit. Needless to say, the argument is identi- 
cal to the one used before for the eigenfunctions living 
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in islands around elliptic orbits, save for = 1 in those 
regular cases. Since orbits of low period are most likely 
to be resolved first, the eigenvalues found for l,nax = 40 
in table || have phases according to p = 1,2 and 4. 

Knowing which orbits are linked to a non-unimodular 
eigenvalue, we can adopt a cycle expansion to calculate 
decay rates from periodic orbits |0. A cycle expan- 
sion of the spectral determinant, i.e. the characteris- 
tic polynomial of the Frobenius-Perron operator, allows 
for the calculation of resonances in hyperbolic system 
with high accuracy [Q. The spectral determinant is ex- 
pressed in terms of the traces of the Frobenius-Perron 
operator Tr as d{z) = 11^=1 exp (-^Tr T^") and 
subsequently expanded as a finite polynomial up to some 
order n^ax- Only the first n^ax traces are required for 
the calculation of this polynomial. The traces Tr 7^" are 
calculated by summing over hyperbolic periodic orbits of 
length n as Tr T'" = ^ | det(i-j)| '^here the 2x2 matrix 
J = dM"'{X)/dX is the linearized map M" evaluated 
at any of the points of a contributing period-n orbit and 
X = {cos9,(p) the phase space point. The zeros of the 
polynomial which are insensitive against an increase of 
nmax are inverses of resonances. 

The condition under which the ordinary cycle ex- 
pansion of a spectral determinant converges is that all 
periodic orbits are hyperbolic and sufficiently unstable 
P, pTp^ . But if we only consider one ergodic region in 
phase space at a time, i.e. bar contributions from el- 
liptic orbits, and impose a stability bound by including 
only the relatively few hyperbolic orbits identified in an 
eigenfunction of T't^), we can still use the cycle expan- 
sion as follows. We assume the spectral determinant fac- 
torized as d{z) = Ili^i ^^ii^) with one factor di for the 
family of eigenfunctions to which a given set of periodic 
orbits contributes. Each such factor di{z) is then calcu- 
lated separately with the above well known expressions 
but restricting the periodic-orbit sum for the trace Tr 
to the orbits previously identified as contributing to the 
eigenfunctions. In table |l| resonances reproduced via the 
spectral determinant from only few orbits show surpris- 
ingly good agreement with the resonance-eigenvalues of 
for r = 10 and Imax — 60. The index rimax gives the 
order up to which the spectral determinant is expanded, 
i.e. the length of the longest (pseudo-) orbits employed. 
The total number of orbits used is given in brackets be- 
hind the resonances. For the resonance 0.8103 the three 
relevant orbits of period 1,2 and 4 are marked in the 
magnified region of the eigenfunction {Imax = 60) in fig- 
ure The first repetition of the single period- 2-orbit 
contributing to the resonance 0.6597 gives an almost di- 
verging contribution to the spectral determinant, thus 
hindering its expansion to a higher order. 

In the cycle expansion the phases of the resonances 
are reproduced exactly since they are again directly de- 
termined by lengths of orbits. If -p is the shortest orbit 



length used in di{z)^ the polynomial can as well be writ- 
ten as a polynomial in thus allowing the zeros to have 
the phases of the pth roots of unity. 

As a final check on the physical meaning of our frozen 
eigenvalues with moduli smaller than unity as Frobenius- 
Perron resonances we have compared these moduli with 
rates of correlation decay. In a numerical experiment 
we investigated the decay of the correlator C{n) — 
[(p(n)p(O)) - (p(cx))p(0))] ■ [(p(O)p(O)) - (p(oo)p(0))]-i. 
Depending on the choice of /9(0) different long-time de- 
cays are observable. We chose p(0) as covering the re- 
gions where the hyperbolic orbits relevant for a given 
resonance are situated. Figure |]c; illustrates the very 
good agreement between the long-time decay of C{n) 
(dots) and the decay as predicted by the corresponding 
resonance 0.81 (full line). Together with the resonances 
at Imax — 60 tabic |l| displays the associated decay fac- 
tors by which C(n) decreases over one timestep, obtained 
from a numerical fit. Again the agreement is convincing. 

In conclusion, we have presented a method to de- 
termine Frobenius-Perron resonances and the associated 
phase-space structures, applicable to systems with mixed 
phase spaces. The acquired knowledge of phase space 
structures allows to check the accuracy to which reso- 
nances are determined by the otherwise independent ap- 
proaches of cycle expansion and correlation decay. 
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FIG. 1. (a) The eigenfunction to the ahnost unimodular eigenvalue 0.9993 for r = 4, Imax ~ 60 is sharply localized on elliptic 
islands surrounding a period-3-orbit in phase space, (b) Phase space portrait of the elliptic islands supporting the eigenfunction. 
(c) The decay of the correlator C(n) (dots) [with the initial density localized in the region shown in figure ^] and the decay 
predicted by the corresponding resonance 0.8103 (full line) agree well (see also table 
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FIG. 2. Eigenfunction to the resonance-eigenvalue 0.81 for r = 10 and (a) Imax ~ 40, (b) Imax ~ 60. While the coarse 
structures are identical, finer structures appear in the eigenfunction with the higher resolution, (c) Magnification of the region 
of large amplitude in figure The numbers indicate the positions of the three periodic orbits of lengths 1, 2 and 4 which are 
used in the cycle expansion up to order Umax ~ 4. 





^max — 50 


^max — 60 


/ — 70 

''max — ' ^ 


0.8116 


0.8205 


0.8103 


0.8076 


0.7457 


0.7547 


0.7470 


0.7459 


-0.7432 


-0.7475 


-0.7510 


-0.7465 


-0.0063 


-0.0123 


-0.0079 


-0.0042 


±i 0.7431 


±10.7515 


±i 0.7517 


±i 0.7414 


-0.6443 


-0.6188 


-0.6377 


-0.6347 



TABLE I. Resonances appearing for r = 10 at Imax ~ 40, 50, 60 and 70. 



^max — 60 


0.8103 


-0.7510 


0.6597 


'^max 1 


0.2185{1} 






"^max — 2 


0.7070{2} 




0.4969{1} 


'^max — 4 


0.7664{3} 


-0.7483{4} 




decay of C{n) 


0.8005 


0.7697 


0.6783 



TABLE II. First row: resonances obtained from the truncated propagator for r = 10 and Imax ~ 60. Below: corresponding 
results from cycle expansion up to order Umax- The total number of primitive orbits employed is given in curly brackets. Last 
row: associated decay factors by which C(n) decreases over one timestep, obtained from numerical fit. 
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